Open-flow mixing: Experimental evidence for strange eigenmodes 
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We investigate experimentally the mixing dynamics of a blob of dye in a channel flow with a finite 
stirring region undergoing chaotic advection. We study the homogenization of dye in two variants 
of an eggbeater stirring protocol that differ in the extent of their mixing region. In the first case, 
the mixing region is separated from the side walls of the channel, while in the second it extends to 
the walls. For the first case, we observe the onset of a permanent concentration pattern that repeats 
over time with decaying intensity. A quantitative analysis of the concentration field of dye confirms 
the convergence to a self-similar pattern, akin to the strange eigenmodes previously observed in 
closed fiows. We model this phenomenon using an idealized map, where an analysis of the mixing 
dynamics explains the convergence to an eigenmode. In contrast, for the second case the presence 
of no-slip walls and separation points on the frontier of the mixing region leads to non-self-similar 
mixing dynamics. 

PACS numbers: 47.52. -fj, 05.45.-a 



I. INTRODUCTION 

Many industrial situations require the mixing of very 
viscous fluids, and involve mechanical stirring in order to 
achieve sufficient homogeneity. We investigate here the 
mixing scenario of a low-diffusivity dye in a free-surface 
open flow that crosses a localized stirring region. Rele- 
vant applications concern continuous throughflow indus- 
trial processes, for example the food processing, glass- 
making, or pulp and paper industries. 

Despite its ubiquity, open-flow mixing has received 
comparatively much less attention than its closed-flow 
counterpart in the applied mathematics and physics lit- 
erature. Efficient stirring flows promote chaotic advec- 
tion [TJ 12] in viscous fluids, that is chaotic Lagrangian 
trajectories of point-like particles. The temporal dynam- 
ics of the concentration fleld of a diffusive dye mixed by 
chaotic advection have been addressed in various exper- 
imental and numerical studies of closed flows [31 IH |SJ 
[i [3 [El 13 Uni CH HI- Time-persistent patterns and 
exponential decay rates of fluctuations observed in ex- 
periments and simulations [H [3 [TOl [El [13] have been 
associated with an eigenmode of the advection-diffusion 
operator dubbed strange eigenmode [5] 1141 1151 1161 I17j . 
Non-hyperbolic structures of the phase space may nev- 
ertheless affect mixing dynamics [18]. In particular, we 
have shown in recent experimental work [T, TW that fixed 
solid walls with no-slip boundary conditions can support 
parabolic separation points leading to algebraic mixing 
dynamics in closed flows ^ 1211 l22l [23] . 

A wide class of open-flow mixing devices successfully 
used in microfluidics rely on the repetition of similar pas- 
sive or active mixing elements — such as grooved pat- 
terns [24], serpentine channels [25l [26], or elementary 
micro-pumps [27] — that induce repeated stretching and 
folding events as fluid particles flow through successive 



mixing elements [HI [551 [55] • However, many large-scale 
processes of heavy industry cannot afford the pressure 
drop induced by repeated geometry changes, which is ac- 
ceptable only for microscales. 

In contrast, we study here a free-surface channel flow 
that crosses a single mixing region where stirrers fold 
and stretch passing fluid particles, which is a more real- 
istic device for large-scale processes. In the following, by 
"open-flow mixers" we refer to devices with a single mix- 
ing region. In such systems chaotic advection necessarily 
affects fluid particles for only a finite time, corresponding 
to their stay in the vicinity of the stirring elements. Suc- 
cessive stages of the mixing of an impurity — a blob of 
dye — are shown in Fig.[l] incoming dye escapes progres- 
sively downstream, and fluid particles experience only a 
finite number of stretching and folding steps. As chaos is 
usually characterized by asymptotic indices such as Lya- 
punov exponents [5] , the description of finite-time chaotic 
advection in open flows requires additional concepts. 

Some tools of dynamical systems are nevertheless rele- 
vant to characterize mixing in open flows. Of particular 
importance is the set of periodic orbits that never escape 
the mixing region, called chaotic saddle [30, 31, 32, 33] 
IM] . This set, together with its stable and unstable man- 
ifolds, are the backbone of fluid transport in and out of 
the mixing region |341 I35j . During their stay inside the 
chaotic region, Lagrangian particles successively shadow 
the trajectories of periodic orbits [3S1 I3Z1 I3H] , and finally 
leave the chaotic region close to the unstable manifolds 
of such orbits. In particular, fluid particles with long 
residence time in the mixing region (Fig. mc-d)) trace 
out the unstable manifold Wu of the chaotic saddle [51] . 
Aspects of mixing in open flows such as residence-time 
distributions, fractal dimensions of the mixing pattern, 
or the dynamics of chemical reactions occurring in the 
flow have been related to properties of the chaotic sad- 
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However, except for an early pioneering study by Soni- 
merer et al. [32] , the evolution of the statistical proper- 
ties of the dye concentration field have never been stud- 
ied quantitatively in open chaotic flows. This contrasts 
with closed flows, where much research effort has been in- 
vested in the study of persistent concentration patterns 
and strange eigenmodes [H [HI [15l [HI [17]. In addition, 
much of the earlier work on chaotic mixing in open flows 
has focused on flows more relevant for geophysics than 
industry (see |39j for a review), such as Von Karman 
alleys in the wake of a cylinder [30] [32] [331 iMj • A note- 
worthy exception concerns mobile point vortices in an 
infinite domain studied by Neufeld and Tel [10], and 
Shariff et al. |41j . which share similarities with moving- 
rod devices that we study here. Our paper follows similar 
lines to the closed-flow studies of dye concentration cited 
above. Rather than characterizing kinematic properties 
of chaotic advection, we concentrate on describing quan- 
titatively the evolution of the concentration fleld of a dye 
tracer over successive stirring periods in an open-flow sys- 
tem, shown in Fig. [l] The question of mixing efficiency 
and its characterization is beyond the scope of this paper 
and left for future communications. Rather, we investi- 
gate thoroughly the evolution of the concentration fleld 
resulting from the initial injection of a tracer. 

The paper is organized as follows. In Section [IT] we in- 
troduce our experimental apparatus, designed to repro- 
duce some key ingredients of open-flow industrial mixers. 
In Section |III| we consider the main qualitative features 
of dye patterns. In Section |IV[ we take on the evolu- 
tion of the concentration fleld for a flrst stirring proto- 
col, where mixing takes place far from the channel side 
walls. After a few stirring periods, we observe the on- 
set of a persistent concentration pattern whose moments 
decay exponentially with time, which we attribute to a 
concentration eigenmode. In Section [V| we introduce a 
simplifled 1-D model derived from the baker's map [42, 
to describe chaotic advection in open flows. The onset of 
the concentration eigenmode and its link with the chaotic 
saddle are investigated in Section [VT] for this ideal model. 
In Section |VII| we report on departure from the eigen- 
mode behavior of the concentration fleld in experiments 
with non-hyperbolic chaotic saddles, a case that occurs 
for example when the chaotic saddle extends to the flxed 
no-slip walls of the channel. We conclude in Section [VIII| 
with some flnal remarks. 



II. EXPERIMENTAL SET-UP AND METHODS 

We have designed an experimental device to study mix- 
ing in open flows with two main criteria: (i) a stirring flow 
that promotes chaotic advection and is simple enough to 
draw generic conclusions; and (ii) a device that is real- 
istic enough to be relevant to applications. To satisfy 
the latter criterion, we require the device to have good 
mixing properties, but still be realizable with standard 
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FIG. 1: Stroboscoped pictures of the dye pattern (separated 
by integer multiples of the stirring period), illustrating the 
successive phases of a typical mixing experiment, (a) A blob 
of dye is advected inside a mixing region where two rods stir 
passing fluid, (b) Stirrers catch a part of the blob, which 
is stretched and folded into thin filaments, while a fraction 
of the blob escapes quickly downstream without experienc- 
ing much mixing, (c) and (d) After a few stirring periods, 
filaments have been stretched enough to reach the difi'usive 
Batchelor scale, and intermediate gray levels appear in the 
pattern. From this moment on, a permanent pattern appears 
to set in: patterns in (c) and (d) resemble each other very 
closely, but with a weaker intensity at longer times. 



mechanical components. 

Viscous cane sugar syrup flows in a long glass channel 
of length 2 m (x direction) and width 36 cm (y direction). 
Sugar syrup flUs the channel up to a height h — 5 cm (z 
direction) . We flx the flow rate at 1 L • min~ by impos- 
ing a constant pressure drop at the inflow of the channel, 
which corresponds to a velocity averaged over a channel 
cross-section U ^ 1mm • s~^. In the absence of stirrers, 
the flow is steady, and independent of x after a few vis- 
cous diffusion lengths Uh^/v ^ 1cm from the channel 
entrance. Analytical calculations and numerical simula- 
tions [13] show that this shallow and free-surface flow has 
a much flatter proflle at all heights than a Poiseuille flow, 
and is instead closer to plug flow save close to the side 
walls. 

Our stirring protocol resembles an egg-beater device 
(Fig. [2|. Two cylindrical stirrers of diameter 45 mm 
are plunged into the fluid and rotated at constant lin- 
ear speed by two motors on intersecting circular paths of 
diameter d = 140 mm. The distance between the centers 
of the two trajectories is 105 mm, allowing a signiflcant 
overlap of the two circles. We can vary the stirring fre- 



quency / to tune the mean number of stirring periods 
spent by a particle inside the mixing region N = d/UT, 
where T — 1/ f. N is a, rough measure of the typical 
non-dimensional timescale of chaotic advection for La- 
grangian particles. TV was varied between 4 and 24 to 
study the dependence on stirring frequency of the evolu- 
tion of the concentration field. 

The direction of rotation of the stirrers with respect to 
the flow direction is an important factor. We can rotate 
the rods so that they assist fluid in passing along the 
channel walls (Figs. l2ja) and (c)), where they tend to 
incorporate fluid inside the mixing region while on the 
upper and central part of their trajectory (see arrows 
in Figs, ^a) and (c)). By analogy with swimming, we 
call this protocol the butterfly (BF) protocol. Or we can 
rotate the rods in the opposite direction, so that fluid 
cannot cross the stirring region along the side walls, but 
is forced into a central funnel (Figs. Wb) and (d)). To 
pursue the swimming analogy, we call this protocol the 
breaststroke (BS) protocol. 

We inject with a syringe a spot of low-diffusivity dye 
(Indian ink diluted in sugar syrup) directly beneath the 
surface, far upstream of the stirring rods. With the stir- 
ring frequencies used here, we did not observe any 3-D 
motion of the dye, which remains just beneath the fluid's 
surface for the duration of the experiments. 

The stirring region is visualized from below using a 
2000 X 2000 Redlake camera of 12-bit depth. We vi- 
sualize the whole width of the channel except for a strip 
on each side of width 3 cm that are hidden by the frame 
bearing the channel. The region of interest (see Fig. fll) 
is lit from above in order to directly relate light intensity 
as measured by the camera to the absorption of light by 
dye, hence to dye concentration. We choose the acqui- 
sition frequency of the camera to be a multiple of the 
stirring frequency. A careful calibration of the response 
of the camera to dye absorption allows us to deduce from 
pixel values the concentration fleld of dye in the whole 
image [H]. 



III. MIXING PATTERNS AND THEIR LINK 
WITH THE CHAOTIC SADDLE 

In this section, we describe briefly some qualitative 
features of mixing patterns in our experiments. In the 
Stokes (viscous) flow regime, the presence of the side 
walls sharply dampens the stirring flow. We may thus 
consider the flow as steady and x-independent outside a 
mixing region. The typical size of the mixing region cor- 
responds to the range of the trajectories of the rods, plus 
a damping length that scales with the channel width. We 
shall now focus on this flow region. 

We obtain completely different mixing patterns by sim- 
ply changing the direction of rotation of the rods, as can 
be seen in Fig. [2J However, some features of the time- 
evolution of the concentration pattern — illustrated in 
Fig. IT] for a butterfly protocol — are common to butter- 




FIG. 2: The butterfly ((a) and (c)) and breaststroke ((b) and 
(d)) versions of the eggbeater protocol impose very different 
stirring flows, hence different mixing patterns, (a), (c) But- 
terfly protocol: the rods push fluid along the side walls. This 
results in a layer of free trajectories along each wall that are 
never caught inside the mixing region. Consequently, after 
a few periods the downstream dye pattern is well separated 
from the walls, (b), (d) Breaststroke protocol: the competi- 
tion between the upstream flow and the flow created by the 
rods results in two separation points on the walls. This time, 
the mixing pattern extends over the whole channel width up 
to the walls. 



fly and breaststroke protocols. When the blob of dye ar- 
rives in the mixing region, some parts of the blob escape 
downstream very rapidly (thick black fllaments in the 
upper part of Fig. flTa)). They keep the same concentra- 
tion level as the original concentration of the blob. Thus, 
our mixer fails to mix such elements efficiently. However, 
dye particles remaining in the mixing region (fllaments 
caught by the rods in Fig.jlja)) are stretched and folded 
by the rods. The dynamics of Lagrangian particles are 
thus chaotic in at least some part of the mixing region, 
leading to an exponential separation of neighboring par- 
ticles at a typical rate A [2j . After several stretching and 
folding events, the width of a filament of dye stretched 
at a rate A inside the chaotic mixing region stabilizes at 
the so-called Batchelor length 
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where the effects of compression and diffusion balance, 
with K the diffusivity. Obviously, stretching is not con- 
stant in the mixing region, but the width of a filament de- 
pends only on its recent history of stretching rate |SJ I45j . 
After a few periods, the dye pattern therefore consists of 



gray dye filaments that have reached the Batchelor scale, 
and of white gaps that result from the repeated incorpo- 
ration of undyed upstream fluid (Fig. fljc)). Most strik- 
ingly, this pattern appears to repeat periodically in time 
(Figs. IT|c) and (d)), and also in space up to distortions 
due to the deviation from plug flow in the channel. 

The persistence of the long-time pattern visible in 
Fig. [l] evokes permanent patterns observed in closed 
flows [U |51 35] that have been dubbed strange eigen- 
modes. For long times, the concentration field converges 
to the strange eigenmode, which is the second slowest 
decaying eigenmode of the advection-diffusion operator 
(the first trivial mode corresponds to a nondecaying uni- 
form concentration in closed fiows). This formalism was 
first presented by Pierrehumbert in a 1994 paper |14j . 
and deeper mathematical foundations were given by sub- 
sequent work [ini m]. The strange eigenmode regime 
corresponds to a global decay of concentration fluctua- 
tions [8j [12] , that arises when a box of width w-q holds a 
large number of fllaments coming from different parts of 
the fluid domain. 

We are not aware of any previous observation of 
strange eigenmodes in open flows. In the remainder of 
this paper, we shall examine whether the presence of a 
permanent pattern, as suggested by the above visual in- 
spection, is confirmed by a quantitative analysis of the 
evolution of the concentration field. We shall see that 
butterfiy protocols indeed give rise to permanent pat- 
terns, whereas the spatial distribution of dye keeps evolv- 
ing in the breaststroke case. We shall also explain the 
origin of such a difference. 

Important differences between mixing dynamics in the 
two protocols may already be inferred from direct visual 
observation of Fig. |2] In the case of the butterfly proto- 
col, after a few periods the mixing pattern has a smooth 
boundary on its upstream side and a cusped boundary 
on its downstream side. In addition, lobes detach at 
each period from the mixing region, as the pattern is 
advected downstream by the main flow and sliced by the 
rods (Fig. [l]). Note that the lobes do not occupy the 
whole width of the channel: there exists a layer of free 
trajectories on both sides of the channel that are never 
caught by the rods and "shield" the mixing region from 
the low-stretching regions at the wall. 

In the case of the breaststroke protocol, filaments are 
smoothly elongated in the upper part of the mixing re- 
gion, and the advected pattern consists of roughly paral- 
lel striations of dye and white filaments that extend over 
the whole width of the channel. One of the striking fea- 
tures of the breaststroke pattern is the presence of two 
separation points on the channel walls (represented as 
small disks in Fig.l2](d)) that stem from the competition 
between the main fiow and the rods. Such separations 
points play a special role, as dye filaments emanating 
from the vicinity of the separation points bear a much 
higher concentration value than in the remainder of the 
pattern (Fig.|2];d)). 

Fixed separation points observed in the breaststroke 



protocol are an obvious, but atypical, example of points 
belonging to the chaotic saddle j47j, which is defined as 
the set of periodic orbits of the flow that return to their 
initial position after a flnite number of stirring periods 
T. Periodic orbits are therefore trapped in the vicinity 
of the rods, unlike most trajectories. Since Lagrangian 
trajectories obey Hamiltonian dynamics in 2-D incom- 
pressible flows, orbits of the chaotic saddle can be either 
hyperbolic, elliptic, or parabolic orbits. Fluid particles 
are stretched exponentially with time in the vicinity of 
hyperbolic points, while they loop around elliptic points 
and experience at most algebraic stretching — elliptic 
points are therefore found in non-chaotic regions, also 
known as elliptic islands. It can be shown [7l [19] that 
separation on fixed no-slip walls are marginally unstable 
parabolic points: fiuid escapes their vicinity with slow, 
algebraic dynamics. Note than the chaotic saddle pro- 
vides an elegant way to define formally the mixing re- 
gion as the smallest connected set that contains all the 
periodic orbits of the saddle. 



IV. BUTTERFLY PROTOCOL: ONSET OF A 
STRANGE EIGENMODE 

In this section, we consider the evolution of the statis- 
tical properties of the concentration field in the butterfiy 
protocol, and we show that the visual observation of a 
concentration eigenmode is confirmed quantitatively. 

We measure the concentration field of dye in a rectan- 
gle inside the mixing region (Fig. pl^a)). The concentra- 
tion mean and standard deviation over the studied do- 
main are plotted for two different stirring frequencies in 
Fig.lstb) and (c). The first thing to note is that the evo- 
lution of the mean is exponential in time, which amounts 
to an exponential distribution of residence times inside 
the mixing region for dye particles injected in the ini- 
tial blob. The exponential evolution of the mean can 
be interpreted as a feature of chaos: a lobe is advected 
away from the mixing region at each period, but chaotic 
advection folds the remaining dye over the entire mix- 
ing region, so that the same fraction of the dye is taken 
away at each period. Residence time distributions with 
exponential tails have been observed in several studies 
of chaotic advection in open fiows [SUl ISl]j at least for 
intermediate times. 

Note that the decay rate is controlled to first order by 
mass conservation across the mixing region. If all fluid 
particles were to undergo a transient chaotic motion and 
be stretched and folded by the rods, the decay rate would 
depend only on the flowrate and the extent of the mixing 
region. In the limit of an infinite stirring frequency, this 
situation would correspond to a "perfect mixer" |48j i4S] , 
in the language of chemical engineering. Perfect mixers 
have the property that a particle leaves the mixing re- 
gion with the same probability distribution as if it were 
picked at random, and have a purely exponential distri- 
bution of residence times. However, our device cannot 
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FIG. 3: (Previous page) Butterfly protocol: (a) Rectangu- 
lar region of the mixing pattern where the concentration field 
is measured, (b) Evolution of the spatial mean and standard 
deviation of the concentration field (in arbitrary units) with 
the number of stirring periods t, for a non-dimensional mean 
stirring time Ai" = 10. Note the similar exponential evolution 
of the first two moments, (c) Same as for (b), for N = 5. 
(d) Histograms of the concentration field at periods 4 to 20 
(A'' = 5). (e) The histograms for periods 7 to 20 collapse onto 
a single curve when the concentration is rescaled by its decay- 
ing mean. Inset in (d): The time cross-correlation of the con- 
centration pattern r(t) = Cov(C(t), C{t + 1)) / {cj{t)ij{t + 1)) 
converges to unity, showing that the spatial pattern repeats 
over time. 



be considered as a perfect mixer as some fluid particles 
escape rapidly downstream along the channel sides, and 
contribute to the short-time peak of the concentration 
mean in Fig. ^b) and (c). The deviation from a purely 
exponential distribution of residence times may be con- 
sidered as a crude measure of the deviation from perfect 
mixing. 

We also observe on Fig. Islb) and (c) that the evo- 
lution of the standard deviation is exponential as well, 
and with the same rate as the mean concentration, a 
first symptom of self-similarity. This is confirmed by 
the plots of the concentration histograms taken at differ- 
ent periods (Fig. p[d)), which collapse on a single curve 



when rescaled by the exponentially-decaying concentra- 
tion mean (Fig. [sFe)). We also observe that the spa- 
tial correlation of patterns taken at successive periods 
(inset in Fig. p[c)) converges rapidly to unity, indicat- 
ing a persistent pattern. This means that after a short 
transient phase the evolution of the concentration field is 
self-similar, and is therefore described by the exponential 
decay of its mean: 

C(x,i) = {C){t) X (7(x) - Coo exp(-t/T)(7(x). 

Such a permanent field is the signature of a coherent pro- 
cess (here, a time-periodic fiow), in contrast to e.g. ran- 
dom turbulent mixing. 

We have thus provided quantitative evidence of a con- 
centration eigenmode in chaotic advection in open fiows. 
Note that an eigenmode scenario was also observed for 
other stirring frequencies for the butterfiy protocol. The 
next section is devoted to an interpretation of the eigen- 
mode. 



V. AN OPEN-FLOW BAKER'S MAP MODEL 

Since Pierrehumbert drew attention to invariant pat- 
terns and strange eigenmodes |14j . much work has been 
devoted to their characterization in closed fiows. The 
mean concentration is conserved in a closed geometry, 
so the evolution of the eigenmode corresponds to an 
exponential decay of fiuctuations around the mean by 
averaging the concentration values of neighboring fil- 
aments [HJ [SD]. It is now widely believed that the 
decay rate arises from a complex combination of spa- 
tial correlations of stretching, and cannot be easily pre- 
dicted [31[ni[Tl[Tl[Tl[5T]. In a recent paper [W, we 
have shown that the unstable manifolds of least-unstable 
periodic orbits are the backbone where the strange eigen- 
mode takes significant values. In the vicinity of these 
orbits, particles repeatedly experience a lower stretching 
than average, and concentration fiuctuations can survive 
longer. 

In this section, we account for the appearance of the 
concentration eigenmode observed in our experiments in 
open fiows. We will rely on a simplified model for an 
open fiow with chaotic advection. Our model is a variant 
of one of the most studied models of chaotic mixing, the 
inhomogeneous area-preserving baker's map fill 1421 1521 

IMI- 

The open-fiow baker's map is defined on a channel- 
like domain, that is an infinite strip (Fig. ffl) composed 
of a square central mixing region, and corresponding up- 
stream and downstream regions. At each iteration, the 
action of the map consists of two parts: (i) a translation 
of the strip by a distance U , mimicking the global ad- 
vection in an open fiow; (ii) a 'traditional' baker's map, 
that is a division into two strips, which are then com- 
pressed and re-stacked as in Fig. [4] preserving area. The 
map has the property of mapping a distribution invari- 
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FIG. 4: (a) The open-flow baker's map: at each iteration, the 
unit square is cut into two horizontal strips, which are then 
stretched and restacked together in a square translated by U . 
Outside the mixing region, fluid is simply translated by U . 
(b) Plot of the one-dimensional open-flow baker's map g for 
7 = 0.6 and U = 0.1. 



ant along the direction y transverse to the channel into 
another such distribution. 

We take our initial "blob of dye" to be a uniform strip 
in the y-direction, i.e. a strip that extends over the whole 
channel width. Inside the mixing region, the baker's map 
stretches and folds this strip to create more strips, leav- 
ing the concentration independent of y |lll I12j . We can 
thus focus on one-dimensional distributions that depend 
only on the x coordinate: they represent a 'cut' across 
a striated pattern such as in Fig. [2|c) or (d). Hence, 
in the following calculation we limit ourselves to a one- 
dimensional version of the baker's map. 

Inside the unit interval, the baker's map g reads 



9 



xi-^ gi{x)Ug2{x), < a: < 1, 



(2a) 



where 



gi{x) = U + jx] g2{x) ^ U + J + {I - -f)x (2b) 



and the union (U) symbol in (2a I means that g is one- 
to-two: every point x has two images given by gi{x) 
and g2{x). The parameter 7 satisfies < 7 < 1 and 
controls the homogeneity of stretching, with 7=1/2 be- 
ing the perfectly homogeneous case. We also define the 
mean stretching rate or Lyapunov exponent, 

logT = 7log7-t-(l -7)log(l-7). 
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FIG. 5: (a) Typical evolution of the concentration profile dur- 
ing the mixing of a blob of dye (left) by the open-flow baker's 
map. Once iterates of the initial strip reach the Batchelor 
scale wb, the pattern repeats over time, with an exponen- 
tially decaying amplitude, (b) The support of the concentra- 
tion profile may be split into two blocks L and R, delineated 
by the iterates of the periodic points that mark the borders of 
the chaotic saddle. The amount of dye inside the two blocks 
obeys simple symbolic dynamics (Eq. Q), which allows to 
obtain easily the decay rate of the profile. 



In the map U represents the flow rate, where the unit of 
time is the period of the map, and the unit of length the 
width of the channel. The baker's map g is represented in 
Fig. |4{b) for the reference values we will use throughout 
the paper, 7 = 0.6 and U = 0.1. 

Under the action of the baker's map, the concentration 
profile evolves as 



C{x,t+l) = C{g-^{x),t). 



(3) 



The map g therefore transforms the concentration profile 
C(x, t) at time t into two images compressed by a factor 
7 and 1 — 7, respectively. In contrast to the classical 
baker's map, a fraction of the strip is nrapped outside 
of the unit interval at each iteration, so that the total 
amount of "dye" in the central interval decreases. Note 
that "leaking" baker's maps relying on similar principles 
have already been used in [40, 54, 55 j however, our map 
follows more closely the geometry of the channel flow. 
Diffusion is mimicked by letting the concentration evolve 
diffusively during a unit time interval |1H I56j . During 
that interval, C evolves according to the heat equation 
with diffusivity k. 



First, we compute numerically the evolution of an ini- 
tial blob under the action of g. The initial condition (left- 
most frame in Fig. IsFa) ) is a strip of constant concentra- 
tion C — 1, located in the interval [0, f7]. A few iterates 
{t = 0,3,7,9) of the concentration profile are shown in 
Fig.lsja) for the parameters 7 = 0.6, U ~ 0.1. The initial 
blob is repeatedly stretched and folded into many thin- 
ner strips (second to fourth frames in Fig. p[a)). Some of 
the strips are mapped downstream, so that the number 
of remaining iterates is less than 2* at time t — or equal 
to 2* for a vanishing flowrate C/ = 0. 

As soon as strips are compressed to the Batchelor 
scale wb, diffusion starts smoothing the strips, so that 
the profiles of neighboring strips are merged together. 
From this moment on, we observe the recurrence of the 
same concentration pattern, but with lower amplitude at 
each iteration (third and fourth plots in Fig.lsla)). This 
eigenmode behavior is confirmed by the same statistical 
analysis as in experiments: as shown in Fig. p[a), we 
measure an exponential decay of the concentration mean 
inside the mixing region, and a similar decay of higher 
moments. Moreover, concentration histograms collapse 
onto a single curve when the concentration is rescaled by 
its mean value (Fig.pTb)). The map therefore reproduces 
well the salient features observed in the physical flow. 



VI. SYMBOLIC DYNAMICS OF THE 
OPEN-FLOW BAKER'S MAP 



We now exploit the simplicity of the map to gain more 
insight into the structure of the asymptotic concentration 
pattern, and to describe its evolution with time. Periodic 
points of the chaotic saddle play a special role in the 
concentration profile, as they form the backbone of the 
profile. Some iterates of dye particles converge rapidly 
towards periodic orbits and remain close to them. The 
support of the asymptotic profile is therefore the set of 
periodic points broadened by diffusion to the minimal 
scale Wb- In other words, the mixing profile is composed 
of a succession of iterates of the interval [0, U] that have 
entered the mixing region at different times: apart from a 
zero-measure set of periodic orbits — the chaotic saddle 
— each point of the interval can be mapped backward 
in time to a point in the upstream region. Iterates of 
[0, U] wider than wb do not overlap with other iterates 
and appear as gaps in the pattern (no periodic points are 
present in these intervals), while iterates that have been 
compressed down to wb overlap with images of the initial 
strip and belong to regions where the concentration is 
positive — the fractal chaotic saddle [33] . 

The central mixing region encloses the support of the 
chaotic saddle, and extends from x — X[ to x — x^, X{ < 
Xc (see Fig. [s] (b)). Here, Xf = U/{1 - 7) = 0.25 is the 
unique fixed point of the map, and Xc is the rightmost 
periodic point in the chaotic saddle. For our parameter 
values, Xc belongs to the unique period-2 orbit of the 
map, Xc = (t/(7 - 2) -H 7)/(l + 7(7 - 1)) ~ 0.974. 



We now look for a partition of the unit interval into 
disjoint regions that are mapped onto each other by the 
map — that is, we search for a Markov partition of the 
map [S71 [5S] , which describes succinctly the dynamics of 
the map. The partition is particularly simple in our ex- 
ample: using the map expression (I2]), one easily checks 
that defining respectively the blocks L and R as the inter- 
vals [x[, gi(xc)] and [g2{xf),Xc\ divides the concentration 
pattern in two disjoint blocks. Moreover, gi compresses 
both L and R inside L, whereas 52 moves a compressed 
image of L to the location of R, and advects R down- 
stream. Defining as Lt and Rt the amount of dye remain- 
ing in the two blocks at time t, we obtain the symbolic 
dynamics [54|l59] 

Lt+i =l[Lt + Rt] , (4a) 

Rt+i = (1 - 7)it. (4b) 

The symbolic dynamics (HI) explains the exponential de- 
cay of the concentration mean: after a short transient, 
the total amount of dye Lt + Rt decays exponentially at 
a rate logr, where 



r = i (7 + v/47 - 372) ~ 0.874 



(5) 



is the largest eigenvalue of the linear system Q. This 
is indeed the rate of decay measured in simulations 
(Fig-ia)). 

Symbolic dynamics therefore allows us to determine 
the exact form of the global asymptotic evolution of the 
concentration profile in the map. Note that the value of 
r depends only on the specific form of symbolic dynam- 
ics ^. In particular, r does not depend on U inside the 
interval where the symbolic dynamics (Hi) is valid. For ex- 
ample, for a fixed value of 7, the symbolic dynamics (Hi 
accounts for the map g in the interval 



(1-7)^^,(1-7)^ 



27 



2-7 



(6) 



We have checked numerically that the dependence of r 
on U follows a devil's staircase, that is a function con- 
stant almost everywhere except at values where the form 
of the symbolic dynamics change. This happens when the 
period of the rightmost periodic orbit Xc changes, as the 
rightmost iterate of the interval [0, U] gradually slides to 
the right as U is increased, until it leaves completely the 
unit interval and it is replaced by another iterate with 
different symbolic dynamics. There is no obvious corre- 
lation in the successive periods of the rightmost orbits as 
U is increased. Symbolic dynamics where the rightmost 
orbit has a short period are valid for a wider range of U 
values than for longer periods. 

The global exponential decay of the total concentra- 
tion value also imposes the same exponential decay ev- 
erywhere in the profile, hence the onset of a strange eigen- 
mode. This can be explained as follows. Choose a box 
centered on x of size given by the local scale of variation 
of the concentration profile wb, and suppose we measure 




7 I \ \ \ L- 

20 40 60 80 
t 



100 



2 4 6 
X = C/{C) 



FIG. 6: Open-flow baker's map: (a) Evolution of the first moments of the concentration field Mi{t) = \J {C(x) — {C)ydx\ 
under the action of the open-flow baiter's map. The mean concentration decays exponentially almost from the injection of the 
dye, while the other moments adopt a similar exponential evolution after a short delay, (b) Concentration histograms at periods 
15, 25, • • • , 95, collapse very well on a single curve. 



a positive value of the concentration C{x) in this box. Af- 
ter a typical number of periods n — logF/log wb (which 
is the number of iterations needed to compress the whole 
unit interval to the Batchelor scale), the concentration 
value is the average of many images of the initial pro- 
file. Neglecting round-off and truncation errors, we may 
consider that C measured at time t represents the to- 
tal amount of dye inside the whole unit interval at time 
t — n, multiplied by the compression factor F", that is a 
quantity 

C{x,t)^r\Lt^n + Rt-n)- 

In the next period t + I, the value of the concentration 
measured at x is now 

C{X, t+l)~ T"{Lt+l-n + Rt+l-n), 

that is 



C(x,t+1) =rC{x,t) 



(7) 



as we have shown that Lt and Rt have a geometric pro- 
gression. Eq. ([7| is exactly the equation of an eigen- 
mode (or invariant distribution) of the advection opera- 
tor (the map). While not rigorous, the above reasoning 
shows that a concentration eigenmode appears because 
of the coherent cascading of large blocks towards smaller 
scales resulting from the repeated compression by the 
same map, whereas the decay of the total concentration 
in the larger blocks is explained by considering symbolic 
dynamics that describe the map. 

We have chosen a set of parameters where the Markov 
partition of the map, and consequent symbolics dynam- 
ics, are very simple. Nevertheless, the same line of rea- 
soning can be followed for more complicated partitions 
to explain how a concentration eigenmodes sets in. In 
the experiments, we know from the study of surface dif- 
feomorphisms [BD] that it is possible to define a Markov 
partition with a finite number of blocks for the real flow 
as well — in the absence of elliptical islands. Hence, we 
can infer the onset of an eigenmode from the existence 
of such a partition. However, it is very difficult to find 



the location of periodic orbits and the Markov partition 
from experimental pictures. In contrast, for the map, 
the partition can be determined by partitioning the unit 
interval using the iterates of periodic points [BT] . 



VII. DEVIATIONS FROM THE STRANGE 
EIGENMODE 

In sections |IV| and |V] we showed the appearance of a 
concentration eigenmode for experiments conducted with 
the butterfly protocol, where the chaotic saddle is located 
far from the channel walls, and we used an open-flow 
baker's map to provide a simple physical picture of this 
evolution. We now report on experiments where the con- 
centration pattern does not converge to an eigenmode, 
which we attribute to the non-hyperbolic character of 
the chaotic saddle. 

We flrst present results for the evolution of the concen- 
tration fleld in breaststroke experiments (see Figs, ^h) 
and (d) for a description of the rod motion). As for 
the butterfly protocol, we study the evolution of the first 
moments and histograms of the concentration, respec- 
tively plotted in Figs. ItFc) and (d) for a typical experi- 
ment. The evolution of the concentration mean plotted 
in Fig. W[c) is close to exponential, though somewhat 
faster. This means that dye increasingly concentrates in 
a zone where it is rapidly advected downstream by the 
fiow. However, the standard deviation, also plotted in 
Fig.lrFc), does not follow the evolution of the mean: it is 
slower than an exponential law. In addition, concentra- 
tion histograms taken at different periods do not collapse 
onto a single curve: Fig. [Wd) reveals that the culprit is 
the high-concentration tail of the histogram, whose rel- 
ative importance grows with time. This is the signature 
of dark filaments reinjected along the unstable manifolds 
of separation points (see Sec. Ill for a discussion of sep- 
aration points) . 

Hence, in the breaststroke protocol the concentration 
is not an eigenmode, since it is neither repeating nor 
exponentially-decaying. This can be related as follows to 



the concentration patterns discussed in Sec. Ill We have 
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FIG. 7: [Color online] Breaststroke protocol: Concentra- 
tion pattern (a) 12 and (b) 15 periods after the initial blob 
enters the mixing region. Beyond the obvious loss of dye, the 
shape of the two patterns differs slightly. The dye approaches 
the unstable manifold of the two separation points and is in- 
creasingly concentrated in the vicinity of the manifold, (c) 
Time-evolution of the mean and standard deviation of the 
concentration field (in arbitrary units). In contrast to the 
butterfly protocol, the moments evolve differently and non- 
exponentially. (d) Concentration histograms, plotted here for 
periods 14 to 34, cannot be rescaled on a single curve as for 
the butterfly protocol. The concentration field therefore does 
not converge to an eigenmode for the duration of an experi- 
ment. 



fluid is slowly reinjected along the unstable manifold of 
separation points, forming dark filaments clearly visible 
in Fig. |7] As time goes on, dye is almost exhausted in- 
side the bulk, whereas there remains a pool of poorly 
mixed fluid near the separation points that slowly leaks 
into the bulk. This contrast between bulk and reinjected 
filaments intensifies with time (Figs. u[a) and (b)), in 
contradiction with the hypothesis of a decaying eigen- 
mode. The discrepancy with a concentration eigenmode 
may therefore be traced back to the differences in the 
stretching rate near the walls, where stretching vanishes, 
and the bulk. 

Structures of the phase space other than parabolic sep- 
aration points may cause departure from an eigenmode 
behavior of the concentration pattern. It is well known 
that elliptical islands can store poorly-mixed fluid in the 
vicinity of their boundary, a phenomenon known as stick- 
yness of islands [S^IISJ. This stickiness is associated with 
exceptionally- low stretching values near the island. For 
the butterfly protocol, increasing the stirring frequency 
causes two small islands to appear inside the circles trace 
out by the rods (Fig. [8|. These stable islands are de- 
stroyed by the main flow at lower stirring frequency. We 
indeed observe that dye accumulates around the islands 
with time, whereas it is depleted in the chaotic bulk. Dye 
is however slowly fed from the islands into the bulk by fll- 
aments escaping from the islands, visible in Fig. [8] This 
evolution of the pattern is incompatible with the onset 
of a strange eigenmode: the evolutions of the mean and 
standard deviation of the concentration are not similar. 



VIII. CONCLUSIONS 



already observed that mixing patterns for the breast- 
stroke protocol differ significantly from those for the but- 
terfly protocol. The patterns extend over the entire chan- 
nel width, and the boundary between the upstream re- 
gion and the mixing pattern is flxed by two separation 
points and their unstable manifolds. Following earlier 
work by Chertkov et al. [20l |21], we have described in 
recent work [TJ [T^ how Lagrangian particles come closer 
to, or escape from, separation points obeying slow alge- 
braic dynamics. As a consequence of the no-slip condi- 
tions, separation points on the walls are parabolic points 
that are only marginally unstable |62j and therefore slow 
down fluid in their vicinity [TJ [12]. Particles approach 
the wall according to slow algebraic dynamics d{t) ^ \/t 
for asymptotic times. A hint of such slow dynamics is ev- 
ident in successive pictures of the mixing pattern, which 
grows slowly towards the unstable manifold of the separa- 
tion points (Figs.lrFa) and (b)). In contrast, for butterfly 
protocols the mixing pattern rapidly adopts a constant 
shape. 

Dye filaments from the initial blob are therefore 
trapped for long times in the vicinity of separation points, 
where they are hardly stretched and mixed. Unmixed 



We have studied the mixing of a dye tracer by chaotic 
advection in an open-fiow geometry, where a net free- 
surface flow crosses a limited stirring region. Our de- 
vice was designed to be relevant for industrial applica- 
tions, but our foremost objective was to classify mixing 
protocols according to generic mixing scenarios. For a 
flrst class of mixers characterized by the absence of reg- 
ular islands in the mixing pattern or separation points 
on the walls, we have observed permanent concentra- 
tion patterns that decay exponentially with time, and 
uniformly in space. A quantitative analysis of the con- 
centration fleld conflrmed that the same pattern repeats 
over time with decaying intensity. The onset of a per- 
manent pattern may be attributed to the convergence of 
the concentration field to an eigenmode of the advection- 
diffusion operator, as it is now commonly accepted for 
closed flows [HI. We have provided here a more physi- 
cal picture of the eigenmode by introducing an open-flow 
variant of the baker's map. Symbolic dynamics of the 
map shed further light on the global decay of the total 
concentration in the mixing region. The repeated and 
coherent compression of the whole concentration proflle 
down to smaller scales then explains the global decay at 
the same rate everywhere in the proflle, hence the onset 
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FIG. 8: (a) Mixing pattern for a high stirring speed relative 
to the mean velocity of the flow. Two elliptical islands are 
visible close to the center of the circular paths of the rods. 
"Branches" of higher contrast emanate from the islands as 
a consequence of their stickyness. (b) The time-evolution of 
the mean and standard deviation of the concentration field 
are exponential, but with different exponents for the two mo- 
ments. This departure from an eigenmode behavior may be 
attributed to the presence of the islands, which retain dye in 
their vicinity for long times. 



of a permanent pattern. 

However, we observed a departure from the eigenmode 
in the case where stretching is anomalously slow in some 
regions of the phase space, that is when either parabolic 
points on no-slip walls or elliptic islands are present in 
the mixing region. The former case is related to a recent 
study [71 [12] where the influence of separation points on 
the dynamics of mixing in closed flows was investigated. 
In the case of a no-slip boundary condition on fixed walls, 
separation points promote slow reinjection of fluid that 
has been poorly stirred in the vicinity of the walls. This 
leads to a net slowdown of mixing dynamics, yielding 



algebraic decay instead of the exponential decay expected 
for a purely hyperbolic system. 

Still, many important questions remain, in particu- 
lar regarding mixing dynamics in the breaststroke case. 
While a previous study [30. reported a power-law shape 
of the residence-time distribution in a Von Karman alley 
in the wake of a flxed cylinder, we have always observed a 
decay of the mean concentration faster than a power law. 
Furthermore, we found the shape of the time-evolution of 
the mean concentration to depend on whether the blob 
of dye had been injected in the center of the channel or 
close to a side wall. Further investigation and modeling 
of this class of mixers is therefore needed. 

A crucial question in chemical engineering concerns the 
quantitative evaluation of the mixing properties of differ- 
ent devices. We have studied here the decay of the first 
moments of the concentration field, which provide infor- 
mation about the time during which a "contaminant" 
remains in the mixing region with a significant concen- 
tration. Other situations of interest for industrial pro- 
cesses include a continuous injection of inhomogeneous 
fiuid (say a succession of blobs of dye) , where one expects 
the output to be "homogeneous enough" on a large range 
of scales. Since the work of Danckwerts [H], a measure 
of mixing commonly used is the intensity of segregation, 
which measures the ability of a mixer to reduce fiuctu- 
ations (as evaluated by the variance of the concentra- 
tion field) between the input and the output. However, 
this measure depends strongly on the spatial distribu- 
tion of the inhomogeneity. In contrast, for the butterfiy 
case, the strange eigenmode characterizes the mixer it- 
self, independent of initial condition. It may therefore 
be tempting to relate the characteristics of the strange 
eigenmode, i.e. its probability distribution function or 
its power spectrum, to the mixer's capacity for reduc- 
ing the concentration fluctuations. Future work will deal 
with the characterization of the early transient pattern, 
which consists of the worst-mixed fluid elements, and ad- 
dress the short-residence-time regime and a subsequent 
assessment of mixing quality. 
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